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Abstract. The correlation calculation of the electronic structure of PbH is 
carried out with the Generalized Relativistic Effective Core Potential (GRECP) and 
MultiReference single- and Double-excitation Configuration Interaction (MRD-CI) 
methods. The 22-electron GRECP for Pb is used and the outer core 5s, 5p and 
5d pseudospinors are frozen using the level-shift technique, so only five external 
electrons of PbH are correlated. A new configuration selection scheme with respect 
to the relativistic multireference states is employed in the framework of the MRD-CI 
method. The [6,4,3,2] correlation spin-orbit basis set is optimized in the coupled cluster 
calculations on the Pb atom using a recently proposed procedure, in which functions 
in the spin-orbital basis set arc generated from calculations of different ionic states of 
the Pb atom and those functions are considered optimal which provide the stationary 
point for some energy functional. Spectroscopic constants for the two lowest-lying 
electronic states of PbH ( 2 n 1 / 2 , 2 n 3 / 2 ) are found to be in good agreement with the 
experimental data. 
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1. Introduction 

The Relativistic Effective Core Potential (RECP) method has been widely used in recent 
years in the correlation calculations of molecules with heavy and super-heavy atoms (see 
e.g. [Q]). This method allows one to treat explicitly only those electrons which play the 
most important role in the making of chemical bonds, while interactions with the core 
electrons and relativistic effects are described with the help of some effective one-electron 
operator (usually semi-local) in the Hamiltonian. Depending on the desired or reachable 
level of accuracy, one can choose an effective Hamiltonian which provides the minimal 
computational effort for the given accuracy. Accordingly, the accuracy of the method 
used for treating correlations and that of the RECP method should be consistent. From 
that point of view, the PbH molecule is of considerable interest as a testing system. 
In calculations on the Pb atom 0, a comparable level of accuracy for the Generalized 
RECP (GRECP) and correlation calculation methods was observed. PbH is a weakly 
bound system so the accuracy of the correlation calculation on the molecule can be 
expected to be close to the accuracy of the correlation calculation on the Pb atom. 

The GRECP method developed in PNPI allows one to attain a level of accuracy 
which is better than 200 cm' 1 for transition energies in correlation calculations of 
different heavy atoms and their ions [0, f|. Our recent calculations || also show that 
the MRD-CI method with an improved selection scheme and properly chosen basis set 
enables one to reach a level of accuracy for molecular calculations involving heavy atoms 
which is comparable to the accuracy of the GRECP calculations for transition energies 
between electronic terms in those atoms. 

It can also be noted that from valence CI calculations and subsequent calculations 
by the Relativistic Coupled Cluster method with Single and Double excitations (RCC- 
SD) with the same number of the correlated electrons, one can easily calculate 
contributions to the total energy from higher-order cluster amplitudes first of all in order 
to account for nondynamic correlation effects in the valence region. That information 
can be used in subsequent calculations by the RCC-SD method with an increased 
number of correlated electrons to give a significant improvement in the accuracy for 
the description of correlation structure of the system of interest @. 

2. Methods 

A detailed description of the methods used in the present calculations can be found 
elsewhere || ||, U . Here we give only brief description and specifics of each method. 

The GRECP operator contains non-local terms with the projectors on outer- 
core pseudospinors, together with the radially-local terms common for other RECPs. 
The addition of non-local terms is connected with a distinction between the effective 
potentials for outer-core and valence electrons. It should be noted that one-electron 
integral calculations with GRECPs are only slightly more time-consuming than with 
ordinary radially-local RECPs, but this time is negligible in comparison with the time 
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needed for two-electron integral calculations. 

In the MRD-CI method || only those Spin (and Space Symmetry) Adapted 
Functions (SAFs) are included in the CI calculation which provide the values of the 
second-order perturbation theory corrections to the energy 

AEi = i < > I 2 (1) 

hi — 

which are greater than some threshold Tq. In the above equation, \l/o is some 
multireference function and are the trial SAFs. As, in particular, was shown 

in || for systems with large Spin-Orbit (SO) interaction, it is preferable to choose \l/o as 
an eigenvector of some Hamiltonian including SO terms. The perturbation corrections 
H for the T=0 threshold to the total state energies are applied after the final SO-CI 
calculations and Davidson corrections ITU are calculated for each term. 



3. Basis set 



To generate and optimize basis set for heavy atoms we mainly followed the procedure 
described in detail in |J. The main features of that procedure are: 

1. A number of the trial one-electron basis functions is generated in Dirac-Fock or 
Hartree-Fock calculations of different ionic states of the atom of interest with the 
frozen core orbitals. In our case, two-component Hartree-Fock calculations with the 
spin-averaged part of the GRECP (Averaged GRECP or AGREP) were carried out. 
As a result, one-electron functions with the same n and / quantum numbers become 
degenerate. Core orbitals are obtained from the AGREP/SCF calculations of the 
Pb 2+ ion (5s,p,d and 6s orbitals in the case of a 22-electron GRECP). Then SCF 
calulations are performed with a different number of electrons in the frozen-core 
orbitals and an electron occupying some valence orbital (e.g. 6p). 

2. The correlation calculation is performed with each generated function included 
in the basis set. We carried out RCC-SD calculations of five lowest-lying 
states of the Pb atom with the 6s 2 6p 2 leading configuration (see for more 
computational details). In these calculations the j-dependent part of the GRECP 
(Effective Generalized Spin-Orbit Potential or EGSOP) is included in the effective 
Hamiltonian together with AGREP. Thus, effects caused by the SO interaction are 
taken into account on the basis set of many-electron spin-orbital functions. It is 
clear that for the adequate description of the correlation effects the number of stored 
nonidentical two-electron integrals on the spin-orbital basis set is substantially 
smaller than that on the spinor basis set. Therefore, calculations on the spin-orbital 
basis sets can be more efficient if correlation is more important than spin-orbit 
interaction as it usually takes place for valence electrons. 

3. The value of some energy functional is calculated using the correlation energies 
obtained in the previous calculations. If the change of this functional (as compared 
to the value of the functional obtained in the previous step) for some function 



GRECP/MRD-CI calculation of PbH 



4 



is maximal in comparison with the cases of using other generated functions and 
exceeds a given threshold, that function is added to the correlation basis set. 
Otherwise, we neglect the function and the procedure starts from step 1 for an 
orbital with other n and / quantum numbers. 



4. Results 



We have carried out three series of MRD-CI calculations for the four states going to 
the two lowest dissociation limits of the PbH molecule. Five external electrons are 
correlated and the 5si/ 2 , ^>Vi/2,z/2 and 5^3/2,5/2 pseudospinors are frozen from the two- 
component GRECP/SCF calculation of the ground state of the Pb 2+ ion (configuration 
6sy 2 ). At the first run, a few of the lowest-energy configuration functions in irreducible 
representations of the nonrelativistic symmetry group were taken for each electronic term 
at each point of the potential curves as the initial approximation for \l/o (i-e. reference 
configurations). For the next two runs, those configurations are chosen as the reference 
configurations which give the largest absolute values of CI coefficients, Cj, at the 
previous run. The number of the reference configurations is selected using "CNFSORT" 
code H such that their total contribution to the wavefunction (C^ = Yli&ef * s 
approximately equal to 96% at the second run and to 98% at the third run for each of 
the considered states and internuclear distances of the PbH molecule (for more details, 
see paper by Mosyagin et al. from the same volume and paper ||). 

As "MRD-CI" program package only allows the use of symmetry group and 
lower orders, all calculations have been performed in the relativistic double group 
symmetry C^. Molecular spin-orbitals are obtained from SCF calculations performed 



by program module "SCF" in "MOLCAS" package QTO | . The correlation spin-orbital 



basis set for Pb [6,4,3,2] generated in the way described above and basis set from 
reduced to [3,2,1] for H are used. We estimate Basis Set Superposition Error (BSSE) by 
calculating the Pb atom in the molecular basis set, i.e. with "ghost" H atom at different 
distances. In our calculations the BSSE does not exceed 70 cm -1 . The following formula 
is used for the T=0 energy correction: 

E T=0 ~ E T=Tk - A AE i > ( 2 ) 

AE,<T fc 

where Ex=T k and Et=q are the total state energies calculated for the Tt and 
zero thresholds, AEj is the energy lowering for the unselected SAFs estimated by 
equation ([I]). The A constant is determined from calculations with two different 
thresholds, T\ and T 2 : 

E T=Tl - A ^] AEi = E T=T2 - A J] AEj . (3) 

A£j<Ti A£ 7 <T 2 

When applying the perturbation energy correction, the value of A never exceeds 3 
with the thresholds from 0.5nH to 0.03/iif. Potential curves for the four lowest-lying 
electronic states are presented on Figure 1. Molecular spectroscopic constants were 
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calculated by the Dunham method in the Born-Oppenheimer approximation using the 
DUNHAM-SPECTR code @. 

Spectroscopic constants for the two lowest-lying molecular electronic terms are given 
in Table 1. The next two electronic levels can not be characterized in the same way, as 
can be seen from Figure 1. 
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Table 1. Spectroscopical constants (R e — bond length, Dq — dissociation energy, 
w e — vibrational constant, T e — transition energy) for the ground and first excited 
states of PbH. 



Calculated Experimental 
Term values values (T^j 
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~ 6900cm" 1 a) 



a ^ Unpublished data by Fink et al. 



The potential curves from GRECP/5e-MRD-CI calculations of PbH. 
Internuclear distance (X-axis) and total valence energy (Y-axis) are in a.u. 




